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Abstract 

We analyse the available data on the polarized asymmetries A\ for proton, 
neutron and deuteron targets. We use a homogeneous and updated set of 
unpolarized structure functions to derive g\ from A\, and we accurately correct 
for the scaling violations in order to obtain gi(x,Q 2 ) with the same Q 2 for 
all x values. The contribution to the Q 2 evolution of a possible large gluon 
polarized density is also considered. The implications for the Ellis-Jaffe and for 
the Bjorken sum rules are discussed. 
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New data on the polarized structure functions g\ of deuterium nuclei® and of 
neutronsil have recently been published, which add to previous data on g\ of protons&S 
Of special importance for the physical interpretation of the results is the derivation 
from the data of the first moments related to the Ellis- Jaffel and the BjorkenB sum 
rules. While the data on gi(x, Q 2 ) are collected at different values of Q 2 for different 
values of x, the evaluation of moments in x requires the same value of Q 2 for all 
values of x. Thus, for a correct evaluation of moments, in principle, one must apply 
corrections to the data points in order to take the Q 2 evolution into account and 
reduce each data point in x to a common Q 2 value. As a first approximation to the 
solution of this problem one can imagine to take advantage of the fact that, within 
the present accuracy of the data, the primary measured quantity, the asymmetry 

Al = an + o-Tt ' W 

shows no appreciable Q 2 dependence 0~S. Then, gx(x, Q 2 ) at fixed Q 2 for all x values is 
obtained from Ai(x) through the relation A\ = gi/F±, with F\ being the unpolarized 
structure function obtained from a fit of existing data. This procedure has been used 
in refs. [0,0. As confirmed by these approximate analyses, the resulting corrections 
to the first moments are small in comparison with the present experimental errors. 
However, in view of more precise forthcoming data, it is interesting to collect the 
known results on the Q 2 evolution of polarized parton densities@0 and describe 
a more accurate method for reducing the set of data to the same Q 2 value for all 
x bins. In fact, an obvious shortcoming of the approximation of refs. is that 
the asymmetry A\ is not predicted to be a constant in Q 2 by the correct evolution 
equations, so that, by assuming it to be a constant, one makes an error of the same 
order as the effect under study. As the whole Q 2 correction is small with respect to 
the present experimental errors, it is no surprise that the variation of A\ with Q 2 is 
not discernible in the present data. 

We describe in the following a more correct procedure that starts from A±(xi, Qf), 
for the 2-th experimental point, constructs gi(xi,Q 2 ) = Ai(xi,Q 2 )Fi(xi,Q 2 ) and fi- 
nally evolves g^x^Qf) into gi(xi,Ql), with Ql being a suitable common value, for 
example the average Q 2 value in the experiment. In passing, we explicitly compute 
the predicted Q 2 dependence of the asymmetry A\ for p, n and d targets, and check 
that the available data are indeed compatible with it. The Q 2 evolution of polarized 
quark densities also depends on the gluon polarized density. We compute the evolu- 
tion for two limiting cases. In the first case we assume a negligible amount of polarized 
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gluons. In the second case we start from a polarized gluon density large enough to 
entirely explain in terms of gluons (with neglible polarized strange sea) the violation 
of the Ellis- Jaffe sum rule observed by the EMC experiment on protons, according to 
the mechanism based on the anomaly proposed in refs. |TT],|T2J and further discussed 
in refs. 0-0. 

We consider the polarized structure functions 



v 1 

g{ = - 2 



4 1 1 A 
-Am + -Ad + -As 



9 n i 



9 
'4 

-Ad 
9 



9 



9 



where 



-Am + -As 
9 9 



q- - q- 



(2) 
(3) 

(4) 



Aq = q + + q 

and q + is the density of the quark q with positive (negative) helicity in a pro- 
ton with positive helicity. We recall that, according to the mechanism discussed in 
refs. p4}4T5|1, the effective polarized quark densities defined by eqs. (0) and (|3]) are 
split into a part conserved by the two-loop Q 2 evolution equations, which is expected 
to be closely related to constituent quarks, and the gluon component, which, due to 
the anomaly, does not decouple from the first moment at large Q 2 . 

We have 



9l = -g? ] + g{ 8) + gf\ 
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g[ S) = - [Au + Ad + As] . 
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The evolution equation for Aq is given by 



Q 2 -^Aq = ?*yp. [Aq ® AP qq + 2Ag ® AP qg ] . 
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Here Ag = g + — g is the gluon density function, and 

4 (l + z 2 \ 



AP qq {z 



The convolution product <g> is defined as usual by 

- 1 dz 



Combining eqs. 



and (pTQ[) , we find 
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or 



2 , 



® AP gg + -Ag ® AP 5ff 



This equation can be approximately integrated by observing that 
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where 2nb = (33 — 2n/)/6 = 9/2. For a small displacement in Q 2 , we can neglect the 
Q 2 dependence of g v { n and Ag on the r.h.s. of eq. (|15D , and obtain 



AP™ + -A# ® AP ( 



0(a?). (17) 



Actually, in the numerical computations, we use the two- loop expression of a s (Q 2 ) in 
terms of A. We take Ajts{ n f = 3) = 383ti 2 g MeV, which corresponds to a s (m 2 ) = 
0.118 ±0.007. With the definitions given above, taking into account the definition of 
the + distribution, one obtains 
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An extreme parametrization of Ag(x) has been proposed in ref. [H] for the EMC 
experiment: 

Ag(x) = Cx-°- 3 (l-x) 7 J (20) 
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with the normalization C such that Jq Ag(x)dx = 5. In the following, for each 
experiment we use the same form for Ag(x), but with a rescaled normalization given 
by 

Ag(x)dx = 5 a ^Mc) (21) 



EXP J 



where Q 2 EMC = 10.7 GeV 2 and Q 2 EXP are the average values of Q 2 for the EMC and for 
the experiment under consideration, respectively. This amount of polarized gluons 
is extreme in the sense that all the EMC deviation from the Ellis- Jaffe sum rule is 
attributed to the effect of gluons, while some amount of polarized strange sea is indeed 
plausible. Also, the above polarized gluon distribution nearly saturates the positivity 



bound from the known unpolarized gluon density especially at large x |14 . 

In the original treatment of the dataS~B, in order to obtain the values of g>i, the 
measured asymmetry A\{x) is multiplied by F\(x, Ql), where Ql is the average Q 2 of 
the experiment. In the present context, obviously we cannot tolerate this neglect of an 
important contribution to the Q 2 dependence. As already stated in the introduction, 
we construct g±(xi,Q 2 ) = Ai(xi, Qf)F\(xi, Qf). At this stage we can implement, for 
all experiments, an updated form for the unpolarized structure functions that takes 
also the recent measurements by the NMC0 into account. For this purpose we use 
the fit to the structure function F2 from ref. |16] and the fit to R = F2/(2xF\) — 1 
from ref. [[17]]. Starting from the resulting values of g%(xi, Q 2 ), we evolve each point up 
to gi(xi, Qq). For the evolution we use a fit to the data points gi(xi, Qf) as an input 
for g\(x) in the r.h.s. of eq. fll?]) and consider both the cases with and without gluons. 
The modified experimental points for gx, and the corresponding ones corrected for the 
evolution to Q^, chosen as the average Q 2 of the experiment, are shown in fig. [[]. In 
the same figure we also show a fit to the data, (solid line), corrected for the evolution 
with Ag = 0. Our fits are given by 

g{(x,Ql) = x°' 2 (l -a;) 3 (0.589 + 1.07a; -1.29s 2 - 0.792a; 3 ) (22) 

g^(x,Q 2 ) = x a2 (l -a;) 6 (-0.670 + 15.0a; -46.0a; 2 + 43.4a; 3 ) (23) 

g™(x,Q 2 ) = x a2 (l -a;) 6 (-0.358 - 1.10a; + 12.7a; 2 -16.5a; 3 ). (24) 

By d we actually mean (p + n)/2, as the corresponding nuclear correction factor has 
been taken into account®. The extrapolation at low x outside the measured range is 
also shown, based on fitting only the last two measured points by Bx 2 (i.e. the same 
a as in the previous fit). Similarly the extrapolation at large x, based on D{1 — a;) 3 
fitted on only the last measured point, is also displayed. 
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In fig. H] we compare the computed Q 2 dependence of the measured asymmetry 
for protons with the EMC data@. For each bin in x the two curves refer to the 
computed Q 2 dependence at the lowest and highest point in x. The Q 2 dependence 
of Ai(x, Q 2 ) is evaluated by using the evolution equations for gi(x, Q 2 ) according to 
the above formalism and the already mentioned (Q 2 - dependent) fit to the unpolarized 
structure functions in order to obtain F\{x, Q 2 ). The results with and without gluons 
are plotted separately We see that the difference associated to the gluons in the 
evolution is small. This is because the evolution equations at xq are only sensitive 
to the parton densities at a; > Xq\ the gluon density, although having a large first 
moment, being concentrated at small x, has little influence on the evolution in the 
measured range of x. We confirm that within the present accuracy of the data no 
evidence for the Q 2 dependence can be seen. But we also see that the effects of the 
Q 2 dependence cannot be ignored in a more precise analysis. In fig. |3] the predicted 
Q 2 dependence of the neutron and deuterium asymmetries is shown. 

The first moments of g\ for p, n and d were obtained by three different methods: 

• (a) We sum the experimental points bin by bin, excluding the first and the 
last bins. We then add to this value the integral of the small and large x 
extrapolations over the remaining x region. This way of treating the end points 
is safer, since the functions g\ are rapidly varying there, especially near x = 1. 

• (b) We integrate our fits of the data over the experimental range excluding the 
first and the last bins. The remaining region is treated as in (a). 

• (c) We simply integrate our fits over the full x range. 

The results corresponding to the different integration methods, with and without 
the corrections for the Q 2 scaling violations, and with and without gluons in the 
evolution equations, are collected in table 1 for p, n, d and p — n, the last being 
the combination of relevance for testing the Bjorken sum rule. For protons, the 
EMC@ and SLAcI data have been combined. It can be seen that the differences 
between moments computed by different methods are well within the errors quoted 
by the collaborations. We prefer procedure (c) of integrating the overall fit over the 
complete range. In fact the sum bin by bin does not take the differences in the errors 
associated with the different bins into proper account, and the extrapolation based 
on only the last points overemphasizes their significance. Nevertheless, we consider 
it important to keep in mind the size of the changes from one method to the other. 
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I 


II 


III 




(a) 


(b) 


(c) 


(a) 


(b) 


(c) 


(a) 


(b) 


(c) 


I(p) 


0.125 


0.122 


0.119 


0.132 


0.126 


0.119 


0.131 


0.127 


0.119 


1(d) 


0.029 


0.029 


0.031 


0.034 


0.034 


0.035 


0.031 


0.031 


0.034 


I(n) 


-0.023 


-0.023 


-0.022 


-0.023 


-0.023 


-0.023 


-0.022 


-0.022 


-0.020 




Bjorken sum rule 


EMC/SMC 


0.192 


0.185 


0.175 


0.196 


0.185 


0.167 


0.200 


0.191 


0.171 


EMC/E142 


0.148 


0.144 


0.141 


0.155 


0.149 


0.141 


0.153 


0.148 


0.140 


SMC/E142 


0.104 


0.104 


0.106 


0.114 


0.114 


0.115 


0.105 


0.105 


0.108 



Table 1: First moments of gi for p, d and n. The group of columns labelled I 
correspond to uncorrected data, those labelled II (III) to data corrected for evolution 
with Ag = (Ag as in eq. (19)). The columns (a), (b) and (c) refer to the different 
integration methods described in the text. 



The results in table 1 are obtained with the central value A = 383 MeV. The effect 
of varying A in the range given above, 267 to 509 MeV, amounts to a relative error 
of the order of 20 to 40% of the computed correction due to Q 2 evolution. However, 
because of the smallness of the evolution correction, the absolute effect on the first 
moments I(p), 1(d) and I(n) is always less than ±0.0012. According to the above 
discussion, for the first moments I(p), 1(d) and I(n) we take the results from columns 
(c) of table 1. The associated errors are taken to be identical to those quoted in the 
experimental papers, as is appropriate, given our indicative purposes. We obtain the 
following results: 

I(p) = [0.119,0.119,0.119] ±0.018 (0.126) (25) 
1(d) = [0.031,0.035,0.034] ±0.027 (0.0245) (26) 
I(n) = [-0.022,-0.023,-0.020] ±0.011 (-0.022), (27) 

where in square brackets we give in a sequence the raw data, the results from the 
evolution with only quarks, and finally those with quarks and gluons. The values 
in round brackets are the raw data values quoted by the experimental collaborations 
(with the error which is also shown). We see that the effect of the Q 2 evolution is 
small with respect to the ambiguities associated with both the fitting procedure and 
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the experimental errors. 

For the Bjorken difference I(p) — I(n), obtained from the three possible combina- 
tions of experiments, we obtain, in the same notation as above: 

EMC/SMC : [0.175,0.167,0.171] ±0.060 (0.200) (28) 

EMC/E142 : [0.141, 0.141, 0.140] ± 0.022 (0.148) (29) 

SMC/E142 : [0.106,0.115,0.108] ±0.058 (0.093). (30) 

Note that the result on the Bjorken sum rule depends on the amount of polarized 
gluons. One should not be surprised by this effect. It arises because the gluons 
contribute to correct the individual data points and the average Q 2 is different for 
each experiment. Also note that the first moments are independent of Q 2 in leading 
order: once the integral has been computed for a given experiment at Q 2 EXP , within 
the present accuracy, it can be combined with the results of other experiments at 
different Q 2 values. As a consequence, for each entry in eqs. (pE|)-(|50]) we cannot 
specify the relevant Q 2 within the range defined by the two experiments. In fact, 
the predicted difference in value for the Bjorken integral is of order (a s /ir) 2 , which is 
beyond the present accuracy. 

In fig. |] we compare the experimental numbers for the Bjorken sum rules, given 
in eqs. fl28"D -(|30"D, with the Q 2 -dependent theoretical prediction, including corrections 
up to order a s /n or (a s /n) 3 calculated in ref. |18[]. The higher- twist correction, as 
computed in ref. |19|] , is negligible in the Q 2 range shown in fig. |j. 

In conclusion, we have studied the effect of scaling violation on the determination 
of the polarized structure function gi for p, n and d targets, and we have proposed 
a procedure to correct for this effect. Our results for the Bjorken sum rule are sum- 
marized in fig. We can see that our results do not differ much from the results 
quoted by the SMC and E142 collaborations. One important aspect of fig. ^ are the 
large error bars on the Q 2 value associated with each determination. We have taken 
these error bars to be the range of the average Q 2 of each pair of experiments. When 
this uncertainty range is taken into proper account, we see that the experimental 
results are consistent with the theoretical prediction, although on the low side. Un- 
fortunately, we also see that the result obtained using the E142 data over lap with a 
region of moderately small Q 2 , where the theoretical prediction in strongly unstable. 
The combination of EMC and SMC data is in a better position from this point of 
view, but in this case the errors are much larger. 
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Figure Captions 

Fig. 1: Effect of the Q 2 evolution correction on the data for the structure function 
gi(x,Q 2 ). Both the raw data and the data evolved to the average Q 2 
of each experiment (with and without the gluon contribution) are shown. 
The solid curve is our fit of the data evolved without gluon. The dashed 
curves represent the extrapolation of the data outside the measured range, 
as described in the text. 

Fig. 2: Computed Q 2 dependence of the proton asymmetry compared with the 
EMC (squares) and E130-E80 (crosses) data. The low (high) curves corre- 
spond to the lower (upper) edge of the x bin. 

Fig. 3: Computed Q 2 dependence of the deuteron and neutron asymmetry, with 
(a) and without (b) the gluon contribution. 

Fig. 4: Theoretical prediction for the Bjorken sum rule as a function of Q 2 . The 
dotted curves are obtained using the 0(a s ) formula, for A = 267, 383 and 
509 MeV. The solid lines include corrections up to 0(a^), for the same 
values of A. The experimental values are also shown. 



